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Abstract 

It is has been long known that the curved space in the presence of gravitation can be described 
as a non-homogeneous anisotropic medium in flat geometry with different constitutive equations. 
In this article, we show that the eigenpolarizations of such medium can be exactly solved, leading 
to a pseudo-isotropic description of curved vacuum with two refractive index eigenvalues having 
opposite signs, which correspond to forward and backward travel in time. We conclude that for a 
rotating universe, time-reversal symmetry is broken. We also demonstrate the applicability of this 
method to Schwarzschild metric and derive exact forms of refractive index. We derive the subtle 
optical anisotropy of space around a spherically symmetric, non-rotating and uncharged blackhole 
in the form of an elegant closed form expression, and show that the refractive index in such a 
pseudo-isotropic system would be a function of coordinates as well as the direction of propagation. 
Corrections arising from such anisotropy in the bending of light are shown and a simplified system 
of equations for ray-tracing in the equivalent medium of Schwarzschild metric is found. 
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I. INTRODUCTION 



It has been shown [1,2] that the propagation of electromagnetic waves in curved space, 
can be described by a mathematically equivalent anisotropic medium in flat geometry. Re- 
cent uses of coordinate transformation [3] and equivalent medium theory [4] have paved the 
way for applications of the theory of general relativity in the domain of artificial anisotropic 
metamaterials. Further applications of transformation media in optical cloacking and in- 
visibility are described in [5], and in a recent extensive review article on the geometry of 
light by Leonhardt and Philbin [6]. Some prior studies are also summarized in [7]. But to 
date and the best knowledge of the authors, no analysis of the governing electromagnetic 
equations for Schwarzschild metric within equivalent medium theory, has been published. 

Here, the optical anisotropy of curved space is demonstrated by means of a rigorous alge- 
braic analysis. We derive the eigenmodes of propagation and conclude that vacuum exhibits 
a property very much similar to pseudo-isotropic media [8] , but with broken symmetry with 
respect to the waves travelling forward and backward in time. A simple pseudo-isotropic 
medium has different refractive indices along all propagation directions, but exhibits no 
birefringence when standard constitutive relations are used [8-10]. This broken symmetry, 
better known as Sagnac effect, arises from the rotation of spacetime, and its implications 
in relativity and on the propagation of light has been investigated and reviewed by Schleich 
et. al. [11]. 

For this purpose, we start by inserting the constitutive relations into the Maxwell's wave 
equations and obtain the governing equation for eigenpolarizations. This is shown to result 
in a modified normal surface equation for the refractive index eigenvalue. Simplification for 
the pseudo-isotropic behavior gives rise to two different refractive indices with opposite signs, 
which are equal in magnitude for a non-rotating spacetime. We are thus led to the conclusion 
that, the speed of light is dependent on the local geometry of the spacetime, but at the same 
time, the curved space of a rotating metric is differently seen by photons, which propagate in 
opposite directions along the time coordinate at velocities shghtly below and above c. The 
difference is easily seen to be removed for non-rotating metrics. In otherwords, the time- 
reversal symmetry of Maxwell's equations breaks down under rotation. This conclusion is 
in contrast to the famous statement by Richard Feynman (1985) [12, p. 98]: "Every particle 
in nature has an amplitude to move backwards in time, and therefore has an anti-particle... 
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Photons look exactly the same in all respects when they travel backwards in time... so they 
are their own anti-particles..." It is thus here shown that Feynman's conjecture does not 
hold for rotating metrics. 

As examples of applicability of our proposed formulation, we consider Nearly Newtonian 
(better known as linearized Schwarzschild metric in isotropic coordinates), Newtonian, and 
Schwarzschild metrics, which are all spherically symmetric and non-rotating. We derive 
exact closed forms for refractive index in these cases. While the Newtonian metric gives a 
non-homogeneous and isotropic refractive index, the latter metric is shown to result in a 
dependence to the radial coordinate as well as the direction of propagation and is explained 
by a non-homogeneous pseudo-isotropic model. We make a comparison to Einstein's results 
in 1911 and 1915 and show that the latter is yet subject to an important correction term. 

In a review paper written in 1907 Einstein mentions that [13,14] ". . .in the discussion 
that follows, we shall therefore assume the complete physical equivalence of a gravitational 
field and a corresponding acceleration of the reference system" . We here conclude according 
to the equivalence principle, that isotropy of velocity of light does not hold for non-inertial 
frames, which makes the space locally anisotropic, and as a result, the refractive index 
should depend on the direction of propagation. It is worth mentioning here that proposals 
of two new satellite missions, namely Gelileo Galilei (GG) [15], and Matter Wave eXplorer 
of Gravity (MWXG) [16], are devoted to precise interferometric tests of the equivalence 
principle with unprecedented accuracies. 

II. THEORY 

An empty curved spacetime may be seen as an equivalent fiat spacetime with a non- 
homogenous and anisotropic hypothetical dielectric filled everywhere [2], the constitutive 
relations of which between electric E and magnetic H fields in SI units take the form [2,5] 

D = £o [£] E + -w X H (1) 

c 

B = /xo [yu] H w X E 

c 

Here, the symmetric dimensionless tensors of relative permittivity e and permeability // are 
given by 

3 



where g''^ and gij are respectively the contravariant and covariant elements of the metric 
tensor of space, with g being the determinant of 4-metric. Also, the gyration vector w is 
defined as 

w = — {goi} (3) 
5*00 

Usually, non- rotating geometries are described with metrics having g^oi = for i = 1, 2, 3 so 
that the gyration vector w vanishes. In fact, it has been shown that w is proportional to 
J X r, where J is the angular momentum of the rotating metric [17,18]. 

We are now able to show that the electromagnetic curved vacuum behaves as a birefringent 
pseudo-isotropic medium, under the sole approximation that the variations of permittivity 
and permeability tensors due to the local geometry of space-time occurs on a length scale 
much greater than the wavelength, such that these two tensors as well as the gyration vector 
w may be taken effectively as constants. As it is shown below, time-harmonic plane-wave 
solutions may be sought in this case. 

We start by plugging (2) into Maxwell's equations, which gives 

k X E = cuB = ujfj,o [fi] H - kow x E (4) 
k X H = -LuD = -lueq [s] E - kow x H 

Here, uj is the angular frequency, k is the wavevector with ko = ui/c being the freespace 
wavenumber. Dividing both sides of (4) by ko, we get 



(nh + w) xE = r]o [fi] H (5) 

(nh + w) X H = -— [£]E 
Vo 

where h — k/|k| is the unit vector along the propagation vector and r^o is the intrinsic 
impedance of vacuum. Furthermore, n — |k|//co is the index of refraction of vacuum; clearly, 
we have k = nkoh. By defining q = h -|- and G = ?7oH, relations (5) may be rearranged 
as 
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nq X E = [//] G 
nq X G = - [s] E 



(6) 



Notice that q is not a unit vector unless gyration vector vanishes, i.e. w = 0. A zero 
gyration vector usually implies no rotation [17,18], as discussed above, however, it should be 
emphasized that not in all systems of coordinates a non-zero w vector denotes true physical 
rotation. For instance, in the Gullstrand-Painlcvc coordinates [19,20], which is connected to 
the Schwarzschild metric through the so-called rain transformation [21], the gyration vector 
does not vanish although the corresponding black hole is static and non-rotating. In fact, 
the time coordinate in GuUstrand-Painleve coordinates is not the exactly the local proper 
time as measured by the black hole, and is obtained through a Lorentz transformation in a 
frame co-moving with a body being accelerated toward the black hole from rest at infinity. 
The equations (6) may be combined to obtain the governing eigenpolarization equation 

as 



CgG = (7) 
>CG = N-'qx{[£]-^qx(-)} + ^ 

£eE = 

^i.= [^]-'qx{[/.]-^qx(-)} + ^ 
Now, if we define the unit vector s — ci/q, then (7) may be further simplified as 



£F = (8) 

£=[/.]-^sx{[e]-^sx(.)} + ^ 

(sj 

Here, F is either of the field vectors E or G, and furthermore, the modified eigenvalue is 
m (s) = nqis) = n |h + ^w| = [nh -|- w[ = [k -|- A;ow[ //cq. Notice that (8) is applicable to 
both of E or G because of (2). 

Since the permittivity and permeability tensors are symmetric, both can be diagonal- 
ized by a proper local orthogonal transformation in 3-space, or local rotation of system of 
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coordinates. However as it is shown later, the final refractive index may be expressed in 
a rotationally invariant form, meaning that such a transformation is not really necessary 
and does not affect the results. In such case, the normal surface equation for the modified 
refractive index eigenvalue m in the local principal system of coordinates would be [8] 



1 1\2/1 1\2 



4+ 4 (9) 



ijizEy rn? J xHySz m? J ^ yi^ix^z m^/ \Hz£x rn"^ ^ 

1 i \ ^ 1 i\ sisi (I i \ ( I i\ sisi 



+ 1 + , V 

The equation (9) is in general a biquadratic equation for the eigenvalue m, and in general 
has two distinct roots for along all directions of s, but the so-called optical axes. We have 
discussed the conditions under which (9) leads to exactly two optical axes corresponding to 
a biaxial medium, or one optical axis corresponding to a uniaxial medium. However, when 
the condition 

f^x /Az (10) 

^x ^z 

holds, which is easily inferred from (2) for the curved empty space, then (9) allows double 
roots along every possible direction s. This situation is referred to as pseudo-isotropic 
medium, in which the medium has no birefringence but is still anisotropic. For a pseudo- 
isotropic medium, the space looks identical along perpendicular directions everywhere, yet 
different along various propagation directions. Hence, all eigenpolarizations are degenerate 
along all propagation directions, and every polarization satisfying k-B = k- D = B- E = 
D ■ H = together with the constitutive relations, would be an eigenpolarization [8-10]. 

After considerable but straightforward algebra, it is possible to show that the solutions 
of (9) for a pseudo-isotropic medium are given by 

^2 ^xl^yl^z ^xf^yf^z (H) 

H^Sl + HySl + ^ySl S ■ [^] ■ S 

Notice that cyclic permutations of a;, and z indices in (11) give rise to the same result 
when (10) holds. For the case of our interest, where [e] — [//], (11) may be rewritten in the 
rotationally invariant form 
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where |-| represents the determinant operation. In the absence of rotation w=0 as well as 
curved geometry, we have m—n and [e] — [p] — 1. Then the plausible refractive indices 
from (12) are given by n — ±1. Non-trivially, the negative solution corresponds to the 
waves travelling backward in time (Appendix A); some earlier works had postulated negative 
refraction [22], which was later shown to be wrong [23]. This situation states in other 
words that in the absence of gravitational field, the electromagnetic field looks the same 
for photons and anti-photons. In contrast to positrons, anti-photons may be regarded as 
photons traveling backward in time [12]. As we shall see below, the symmetry is broken 
in the presence of rotating gravitational field, and as a result, photons and anti-photons 
become different at least with regard to their velocities. 

The actual solutions of (12) for the eigenvalue n become complicated when we take the 
dependency of m on s into account. By definition, (s) = \nh + w|^ and s — q/q. Then, 
we have 



which by noting q = h -|- can be further simplified to obtain the master eigenvalue 
equation for the refractive index of curved vacuum as 



nh -I- w 




(13) 



(nh + w) ■ [^] ■ (nh w) = 1 



(14) 



Here, [^] — [p] / Expanding (14) in its components gives 



(15) 



which can be rearranged as 



+ 2Bn - C = 



(16) 



S = h • [C] • w 
C = 1 - w • [,^] • w 
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Therefore, the exact solutions are 



VACTW-B \/(h.l5].h)(l-w.l5).w) + (h.K).wf -h-iei-w 

= + A = iTW^ (1^) 



n2 



^/AC + W+B ^(h.[e]-h)(l-w[e]-w) + (h-[e]-wf + h-[e]-w 



A h-[e]-h 

If deviation in metric from the fiat Minkowskian geometry is not too strong, then the 
terms being second order in w can be dropped and (17) can be approximated as 



1 h • [f 1 • w 



h-[e]-h 



Note that (17) take on fairly simple forms when w=0, and thus we get the exact expressions 



ni= + (h.[e]-h)-^ (19) 
n2 = -(h-[e]-h)-^ 

This shows that for non-rotating metrics when described in proper coordinate system, time- 
reversal symmetry holds. It should be also added that a prior study [24] has shown that 
vanishing birefringence is necessary for consistency, where by birefringence the authors of [24] 
mean the existence of two refractive indices for forward-propagation in time, corresponding 
to two light cones. Our study also justifies this result, and furthermore we notice that two 
inequal roots are given by (17), although they may differ only within a sign under some 
circumstances. Hence, there is one light cone for forward propagation and another light 
cone for backward propagation in time. 

Time-reversal of Maxwell's equations plus two orthogonal polarizations allow for four 
distinct refractive index eigenvalues. Notice that for standard anisotropic media, from the 
solution of normal-surface equation we normally have two pairs of eigenvalues, rii — —n^ 
and n2 — —n^ with either ni > n2 > or ^2 > ni > 0, where the equality sign holds only 
along optical axes; it is clear that the time-reversal symmetry holds. This is while in the 
pseudo-isotropic medium of our interest, we have two double roots given by rii — and 

8 



712 = n-4 with either 712 < — 1 < < ni < 1, or — 1 < ^2 < < 1 < ni depending on 
the sense of rotation; hence, in general, the time-reversal symmetry does not hold. This is 
while for non-rotating metrics we have the further simplification rii = —n2, corresponding 
to vanishing birefringence and preserving time-reversal symmetry. Based on the definition, 
a pseudoisotropic medium has no birefringence at all, but the refractive index can still be 
dependent on the direction of propagation. 

When w 7^ 0, we may still note that h is a unit vector, and thus (18) can be still simphfied 
further in the weak gravitational field limit as 

ni^+l + I- An (20) 

712 ~ —1 -|- 7^^^ 

corresponding respectively to photons and anti-photons; here, we define Ari = — 2h ■ [^] ■ w. 
Firstly, it can be seen that the curved vacuum exhibits a local timc-rcvcrsal asymmetry given 
by |A7z| = \n\ — 712], which is roughly a linear function of gravitational potential (as shown 
below) . Secondly, photons travel at a speed shghtly below (above) the speed of light in flat 
vacuum c, while anti-photons travel at a speed slightly above (below) c, if the direction of 
propagation is anti-parallel (along) to the rotation of the universe, where ISn > (A7i < 0). 

Now as it is shown below for Nearly Newtonian, Newtonian, and Schwarzschild metrics, 
the equations (19) are compared to the predictions made by Einstein for weak gravitational 
fields [25-27]. 

III. EXAMPLES 

In this section we consider two important cases, which both correspond to spherically sym- 
metric non-rotating universes: Newtonian and Schwarzschild metrics, where two versions of 

Newtonian metric are discussed. Rotation can be exactly implemented through (17), how- 
ever, the resulting expressions are too complicated and hence are not discussed for the sake 
of convenience. 
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A. Nearly Newtonian and Newtonian Metrics 



The spherically symmetric Nearly Newtonian metric [28, p. 445], or the so-called lin- 
earized Schwarzschild metric in isotropic coordinates is given by the line element 

ds^ = -c^ (l - ^) df + (l + ^) {dx^ + dy' + dz") (21) 

where Vg — 2GM/c^ is referred to as the Schwarzschild radius of the star with M and G 
respectively being its mass and gravitational constant, and dl"^ — dx^ -\- dy^ -\- dz^ is the 
spacelike path element. This metric is an approximate solution of Einstein field equations, 
but to high precision most stars are static and spherical [28, p. 446], so that (21) is applicable. 
Landau and Lifshitz [29] report the Newtonian metric in the slightly different form 

ds'' = -c^ (l - 2^) dt" + {dx" + dy^ + dz") (22) 
The solution (19) may be used to all three above cases and we readily obtain 



r2 (r - Ts) 



1 i 

f2 ( f \ 1 



n = ± Y = ± I 1 - 2 



(r -2r,)2 

respectively for the metrics (21) and (22). Denoting $ = —fs/r, we get 



1 + $ ^ ^ 

1 

n = ±- 



^1 + 2$)2 

These expression blow up to infinity for r — as (r — rs)~^, but approach in magnitude to 
unity as r ^ oo. In the limit of small normalized gravitational potential r >> Ts, however, 
we get 



\n\ 



\n\ 



2 

^ 1 - $ 



(25) 
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The latter result agrees to that of Einstein's 1911 early prediction [25]. A few years later, 
however, he showed that the change in refractive index should have been given by \n\ ~ 1—2$ 
[26,27]. Now, it is discussed below that his correction factor of 2 was still inaccurate. 

B. Schwarzschild Metric 

Schwarzschild metric is known by Birkhoff's 1932 theorem [28, p. 843], to be the 
most general solution of Einstein field equations under spherical symmetry and no rota- 
tion. Schwarzschild metric has been written in many coordinate systems which are all 
connected through transformation, including isotropic, synchronous (Lemaitre-Rylov's), 
Eddington-Finkelstein, Kruskal-Szekeres', harmonic, as well as GuUstrand-Painleve coor- 
dinates [19,20,30]. In its basic form, it is given by [28, p. 607] 



where (r, 6, 0) constitute the standard spherical polar coordinates. Recently, Virbharda has 
revisited the problem of light deflection in Schwarzschild geometry [31-33]. 

The difficulty in working with this metric arises from the fact that the spacelike path 
element df — dr'^ + r'^dfl'^ where dQ'^ — dO"^ + s\v? Odcf)^, does not appear exphcitly in the 
metric. To overcome this difficulty, the common practice is to express (26) in a similar form 
to (21) using the so-called isotropic coordinates [28, p. 840], through the transformation 



where p is called the isotropic radial coordinate. This allows to rewrite the metric as [34] 




(26) 



r 




(27) 



ds^^-c^^ — -^de+[i + ^] [dp'' + p'{de'' + sin^ed^^)] (28) 





which allows one to readily obtain [34,35] 




(29) 



This result can be approximated in the limit of large p as 
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1 + 2— 1 + 2— , p»rs (30) 
p r 

By comparing (30) with (24) the Einstein's correction factor of 2 becomes evident. There 
is no such apparent anisotropy in (30), as obtained from transformation to isotropic coordi- 
nates. 

But the new radial coordinate p in isotropic coordinates has no direct physical meaning, 
while the radial coordinate in (26) is actually the circumferential radius [36,37] (proper 
circumference divided by 27r) . Therefore, we stick to the Schwarzschild coordinates to avoid 
wrong conclusions. However the conversion process to quasi-Minkowskian coordinates needs 
some care, and construction of the proper metric needs the following transformations [36, 
p. 181]. First, the metric (26) is rewritten as 

ds'' = -c^ f 1 - ^) dt" + -^dr"" + df (31) 
V r / r — Tg 

and then we note that r = a/x^ + + z^. This enables us to write down 

_ ^dx'dx^ (32) 
in which ^,j=l,2,3 denote x, y, and z coordinates. Hence (31) can be rewritten as 

ds'' = -c' (l - ^) de + {^-^^ + 8,^ dx'dx^ (33) 

This metric is obviously non-diagonal, but can be diagonalized everwhere through orthogonal 
transformations in 3-space. Since eigenvalues of the resulting tensors are degenerate, Gram- 
Schmidt ortho-normalization was applied to find the principal axes of coordinates. It should 
be noted that the propagation vector h should be rotated accordingly to obtain correct 
results. Furthermore, (19) are exact as w = 0. 

Since h is a unit vector it can be described by the spherical polar angles {tp, x) iii the 
original non-rotated reference frame. A tedious but straightforward calculation then gives 
the eigenvalues of propagation as 




- 1^ - 1 + [cos (0 - x) cos (^ - V') + 2 sin^ (^) cos 9 cos V'] ^ 
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Spherical symmetry makes the absolute choice of the angles (■0, x) irrelevant, in the sense 
that we may set the z-axis along the position coordinate r and choose = ^ = 0. Hence, 
we get the exact yet fairly simple expression 



\n (r) I = ^ (35) 

where the azimuthal angle of propagation if) is measured with respect to the position vector 

_3 _2 

r. Again for r — the refractive index blows up respectively as (r — r^) 2 and (r — r^) , 
if -0 7^ ±1 and ip = However, it approaches in magnitude to unity as r — > oo, since the 
metric relaxes to that of the Minkowskian in the limit of infinite radius. 
Finally, we adopt $ = —Ts/t, and rewrite (35) as 

|n(r)| = (1 + $)"2 (1 + $sin^^/;) 2 (36) 
This result can be rewritten in the limit of small $ (weak gravitational field) as 

l^(r)l~l (37) 

Since is actually the angle made by the position vector r and wavevector k, the above 
equation can be also put in the more convenient form 



|n(r,k) 



fk-r) 



2 



$ = 1 - (2 - ^)$ (38) 



2A;V2 

where hr = kr/k = h ■ r/r. The correction factor to the Einstein's 1915 result as shown in 
(30), hence actually varies between 3/2 and 2 depending on the angle of propagation. The 
light ray evidently keeps passing on the plane defined by the vectors r and k. Then hr — 
or sin^ t/j — 1 holds only at the nearest point to the center of the star in the hght trajectory, 
where (30) is accurate, while at farther points away from the center of the star, we approach 
hr = 1 or sin^ -0 = 0. 

Another conclusion is that this anisotropy is expected to be present everywhere around 
a massive object, so that the change in refractive index by changing the direction of prop- 
agation from = I to = 0, could reach as high as |$| /2 = (GM/c^) /r. Based on the 
estimates given in [25, p. 459], this figure should be of the order of 10~^ for an experi- 
ment done at Earth's distance from Sun, while it would be only about 6 x 10~^° at the 
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surface of Earth when the gravity of Sun is neglected. Finally, the equivalence principle 
states that gravity should be indistinguishable from pure acceleration. Therefore, acceler- 
ated (non-inertial) frames should have optical anisotropy, meaning that the speed of light 
along various directions is not necessarily the same for a non-inertial observer. 

1. Deflection of Light 

Here, we investigate the effect of anisotropy in deflection of light around a massive 
Schwarzschild object. When the coordinates are expressed in the isotropic system of co- 
ordinates, the angle of deflection is simply given by [38-40] 



where ro is the closest distance of approach of the light ray to the star. However, when 
the effects of anisotropy are included, one must solve the associated equations of motion 
in the Schwarzschild metric. A classic solution to this problem in Schwarzschild geometry 
has been proposed by Weinberg [36], and more recently been revisited by Virbhadra [31-33]. 
This results in the expression 



where by comparing to (39), it becomes possible that the first correction term may be 
attributed to the anisotropy of Schwarzschild metric. 

2. Ray Tracing Equations 

Propagation of light rays within the equivalent medium framework can be approximated 
by geometrical optics. This results in two coupled sets of equations, together describing 
the equations of motion for the beam [4,41,42]. The Quasi-Isotropic Approximation (QIA) 
[43,44] for weakly anisotropic media is best suited to deal with this problem when the 
rotation of polarization is sought. This is particularly appealing, since the equivalent medium 
of vacuum is both pseudo-isotropic with vanishing birefringence and weakly anisotropic. 
Anyway, the exact equations of motion take the form 



2r, 



(39) 



a = — 




(40) 
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dr 

d^ 
dk 

d^ 



dH 

dH 
dr 



(41) 



Here, the Hamiltonian H is expressed by if = /(r)(k • [^] • k — |^|), with /(r) being an 
arbitrary function of position. The above set of equations result in six nonhnearly coupled 
first-order differential equations, which normally admit only numerical solution. A great 
deal of simplification is possible for the stationary case when the free parameter r = / is the 
same as arc length covered by the ray [37] 



dr 

dl 



b(r,k) 
|brr,k)| 



(42) 



dk 

H 





k 


dn-\r, k) 




b( 


r, k) 1 dr 



where 



b(r,k) 



d\k\ n-^(r,k) 
dk 



(43) 



Now from (38), we may insert n~^{r, k) 1 + (2 — i/i^)$, which by defining the unit radial 
vector r = r/r results in b(r, k) = h + $ (2 + ^h"^)] — K^f. We also have b{r,k) = 
^b(r, k) • b(r, k), which after removing second- and higher-order terms in $ and applying 
binomial expansion takes the form b{r, k) 1 + $(2 — |/i^). Furthermore, the gradient in 
the right-hand-side of (42) simphfies as 



dn ^(r, k) 
dr 



r 



-/i^ ) f + hrh 



(44) 



Therefore, we obtain the ray tracing equations for the equivalent medium of Schwarzschild 
geometry, being correct to 0($) as 



dk 

'dl 



dx 

dl ' 

-w 

r 



(1 + h - hr^r 
2 



(45) 
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IV. CONCLUSIONS 



In summary, we have discussed the optical anisotropy of vacuum subject to a gravitational 
field. Exact expressions for the refractive indices are obtained, and shown to correspond 
to waves traveling forward and backward in time. These two eigenvalues are different in 
magnitude for a rotating spacetime, meaning that photons and anti-photons look different 
in rotating curved spaces, travelling at velocities slightly below and above c. We also derived 
exact and closed-form expressions of the refractive index for Nearly Newtonian, Newtonian, 
and Schwarzschild metrics, and showed that the latter metric gives rise to a locally pseudo- 
isotropic universe. It has been furthermore proved that a strong local anisotropy is predicted 
to exist in the equivalent medium description of curved geometry, which shows that the 
Einstein's 1915 correction factor might have still been inaccurate. Based on the equivalence 
principle, gravity is indistinguishable from pure acceleration. Hence, accelerated non-inertial 
frames induce optical anisotropy in the speed of light. 
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APPENDIX A: NEGATIVE REFRACTIVE INDEX DOES NOT IMPLY NEGA- 
TIVE REFRACTION 

We here show that existence of a negative refractive index does not imply negative refrac- 
tion. This fact is more than trivial, but can be easily observed by checking the Maxwell's 
equations for simple classical vacuum with flat geometry: 



V X E 



V X H 




(Al) 
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Substitution of solutions varying as exp \j {cut — k • r)] gives 



k X E = +//o<^H 
k X H = —Eqcu'E 

The wave equation for fields may thus be written as 

k X (k X E) = -cu'^iiqEqE 
k X (k X H) = -cu'^f^oCoH 

Noting that from Maxwell's divergence equations we have k • E = k • H = 0, the above may 
be rewritten as 

(k • k) E = uj^HoSoE = — E 
(k • k) H = u;Vo£oH = 
Both equations are satisfied if the following dispersion equation holds 

k • k = A;' = ^ (A2) 

For a fixed direction of k, the two solutions to the above second-order equation are ki = 
ni^s and k2 = n2^s, where s is a unit vector, = = 1, and |-| denotes the absolute 
value. Clearly, we have the acceptable solutions as ki = and k2 = —'^s. Because of the 
obtained forms, we may choose the sign of cu arbitrarily as positive, but will keep it fixed 
throughout. We can show that both solutions do have physical interpretation. 

Note that both solutions to the dispersion equation (A2) are acceptable while there is 
evidently no negative refraction. Now we insert the solutions in the original Maxwell's 
equations to get 

+/xowH ^ s X E = +G (A3) 
-eocoE ^ s X G = -E 
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H — s X E = 

c 

H — s X H = 

c 



s X E = +/io<^H =^ s X E = -G (A4) 

c 

s X H = —Equ'E =^ s X G = +E 

c 

where 770 = c/Iq and G = rjoH. 

Notice that the first set of equations (A3) corresponds to a right-handed coordinate system 
constructed of {s,e,g), while (A4) corresponds to a left-handed coordinate system {s,e,g), 
in which e = E/ |E| and g = G/ |G| = H/ |H| are field unit vectors. This is while the 
left-handed solution (A4) could be obtained equivalently from the ansatz exp [j {cut -|- k • r)], 
too. 

In summary, (Al) admits two different particular solutions, given by 



E (r, t) = Re [Eq exp \j {ut-k- r)]] = |Eo| cos [ujt-'k-r + ^Eq] 
H (r, t) = Re [Hq exp [j {oot-k- r)]] = |Ho| cos [ujt -'k-r + IHq] 

with either Hq = x Eo/770 = k x Eo/a;//o and k = +u;s/c, or Hq = — s x Eo/770 = 
k X 'Eq/uho and k = —tos/c. In both solutions, we may set |Ho| = |Eo| /t^o = 1 and 
ZE'o = Zi^o = for the sake of simplicity. Now, by setting (s, e, g) = {x, y, z) we get 



E (r, t) = y7]o cos uj \^ x 

H(r, t) = icos uj ^ X 

as the right-handed solution. Similarly, we may take (s, e, g) — {x, z, y) to get 



(A5) 



E (r, t) — zrjo cos 
H (r, t) — y cos 



cu { t + -X 
u { t -\ — X 



(A6) 



as the left-handed solution. 

It is easy to verify by direct substitution in (Al), that (A5) and (A6) are indeed both 
correct solutions of Maxwell's equations. By inspection, it may be easily seen that (A5) 
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represents a wave that propagates along +x {—x) when time increases (decreases). At the 
same time, (A6) represents a wave that propagates along +x {—x) when times decreases 
(increases). Hence these two solutions may be better coined as waves traveling forward and 
backward in time, corresponding to the trivial refractive indices Ui = +1 and n2 = — 1, 
respectively, represented as 



El (r, t) = yr]Q cos 







ni \] 




— x 






c / - 






ni y 


UJ 


X 






c /J 



(AT) 



and 



E2 (r, t) = zrjo cos 
H2 (r, t) =y cos 



u (t 



—x] 
c / 

—x] 
c / 



(A8) 



In conclusion, we have the following: 

1) The apparently negative refractive index —1 of simple vacuum, corresponds to the waves 
traveling backward in time. 

2) Both solutions appear as right-handed if the triad is taken as (j;, e, instead of (s, e, g). 

3) Since the unit vector along Poynting's vector lies along p — e x g, then p • A; = 1, or 
P • k > 0. In other words, there is no such "negative refraction" . Notice that P • s > may 
not hold for (A8). 
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